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Abstract. We review the study the Wandzura-Wilczek relation for the structure function g2, with 
a particular attention on the connection with the framework of Transverse Momentum Dependent 
factorization. We emphasize that the relation is broken by two distinct twist-3 terms. In the light of 
these findings, we clarify what can be deduced from the available experimental data on g2, which 
indicate a breaking of the order 2Q-A0%, and how to individually measure the twist-3 terms. 
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At large virtuality Q^, the lepton-nucleon Deep Inelastic Scattering (DIS) cross sec- 
tion scales with x = / (2Mv) modulo logarithmic corrections, where M is the target's 
mass and v the virtual photon energy. At low Q^, power suppressed contributions be- 
come important, e.g., target mass corrections of 0{M^ /Q^) and jet mass corrections of 



order 0{mj/Q ) with nij the invariant mass of the current jet |!lil2|], and higher-twist 



jet ma 



(HT) corrections of C>(Ag(^^/ Q^) related to quark-gluon correlations inside the nucleon 
[jSl]. There are many reasons why we need to identify and measure higher-twist terms 
in experimental data, for example, (i) to verify quark-hadron duality [4], (ii) to measure 
twist-2 Parton Distribution Functions (PDF) at large fractional momentum x and low- 
Q^, e.g., the dju and Ad/d quark ratios, sensitive to the nonperturbative structure of the 
nucleon dl, 0, 13], (iii) to measure multiparton correlations, important to understand the 
nucleon structure beyond the PDFsflsj], (iv) to determine the perturbative QCD evolution 
of the g2 polarized structure function among others [9], (v) to calculate the high-^r tails 
of transverse momentum dependent (TMD) parton distribution functions [10]. 

The inclusive DIS cross section is determined by the hadronic tensor W^^, defined 
as the imaginary part of the forward virtual photon Compton scattering amplitude. W^^ 
can be decomposed in 2 unpolarized structure functions, Fi 2, which we do not discuss 
here, and 2 polarized structure functions gi^2- Its antisymmetric part reads 



F ■ q L V F ■ q J 



(1) 



where F,S are the target momentum and spin. Among the structure functions, g2 is 
unique because it is the only one with twist-3 contributions that can be measured in 
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inclusive DIS. Furthermore, its higher-twist contribution can be isolated thanks to the 
Wandzura-Wilczek (WW) relation, first obtained in the Operator Product Expansion 
(OPE) formalism 



(2) 



Here is determined by the leading twist (LT) part of ^i, which is rather well known 
experimentally: 



Jx y 



(3) 



Strictly speaking, the WW relation is the LT part of ©■ The breaking term A is a pure 
HT term, meaning that its moments are matrix elements of local operators of twist-3 or 
higher, with "twist" defined as dimension minus spin of the local operator lll2l ll3ll. In 
this talk we will limit our analysis to twist-3 operators, and will drop the dependence on 
2^ for ease of notation. Details can be found in Ref. [11411 . 



PARTON DISTRIBUTIONS IN PERTURBATIVE QCD 

In perturbative QCD the structure functions can be expressed as a convolution of a per- 
turbatively calculable coefficient, and a number of nonperturbative LT parton distribu- 
tions and HT parton correlations. In particular, the WW relation can be obtained in the 
framework of coUinear factorization [9] or in transverse momentum dependent factor- 
ization, as we will shortly describe. 
Let's define the quark-quark correlator 



%{xrkT) = I "^iff ^'"'^ S I v>J(0) (0, ^ In-) v/f ) I 5) ^, „ (4) 



{2n 



where /, j are Dirac indices, a is the quark flavor index, k its 4-momentum, x = k ■ 
n-/P-n- its fractional momentum and kj its transverse momentum relative to the parent 
nucleon. The "plus" and "minus" components of a 4-vector are defined as = a • 
in terms of two orthogonal light-cone vectors n\ = = such that n- ■ n+ = 1, and 
n'^ is proportional to up to mass corrections. V7 is a Wilson line (gauge link) whose 
precise form depends on the process. The direction of the Wilson line is determined by 
an additional 4-vector beside P, 5, which in tree-level analyses such as we pursue here is 
identified with the light-cone vector n-. In the light-cone gauge n- -A — the Wilson 
line is identically equal to 1 and 



Nonetheless, the dependence on n_ appears explicitly in the gauge field propagators and 
cannot be in general neglected. 



For any Dirac matrix T we define tlie projection c^"^^^ = Tr[r4>"]/2. The relevant 
TMDs are defined as follows: 

The inclusive DIS is determined by coUinear parton distribution functions (PDFs) which 
are defined by transverse momentum integration of the TMDs: g^{x) = J d^kj g^ix.kr) 

and g[^\x) = f d^kr ^g^{x,kT), with ^ indicating any of the above defined TMDs. 



EQUATIONS OF MOTION AND LORENTZ INVARIANCE 



The Dirac equations of motions for the quarks, and the Lorentz invariance of the theory 
imply the following 2 relations between twist-2 and pure twist-3 functions: 



(EOM) 
(LIR) 



gf^'\x) =xg"Ax) -xg"T{x)+0{m/M) 
m + j-/^\x)+g"r{x) 



sUx) 



(6) 
(7) 



where for light quarks we can neglect the term proportional to the quark mass m 
compared to a typical hadronic scale M. g and g are pure twist-3 parton correlation 
functions (PCF) defined in terms of the quark-gluon-quark correlator, which in the light- 
cone gauge reads 



^jix. x') / .'^-^ {P\ V/, (0) ig d^A-iv) Vii^) \P) 



where a is a transverse index, x' 



^+=^7-=o ' (8) 

J7 +=777-=0 



p and F is the QCD field strength tensor. The 
Lorentz decomposition of defines the relevant PCFs [T5l.[l6ll. 



M 



Gf{x,x')i£p"STa + GF{x,x)S^r5 + ■■■ 



(9) 



where hermiticity and parity constrain Gf{x,x') = Gf{x',x) and Gf{x,x') = —Gf{x',x). 
The pure twist-3 functions in Eqs. ©-(IT]) are particular projections over x' of Gf{x,x') 
and G/7(x,y)(PV denotes the principal value): 



xg'^{x)=VW Jdx 



gT[X) 



PV / dx' 



, Gf{x,x') +Gf{x,x') 

2(y-jc) ' 

, Gf{x,x')/{x — x') 
x — x 



(10) 
(11) 



and as such are sensitive to different parts of the quark-gluon-quark correlator. It is 
very important to find several such quantities, because physically it is only possible 



to measure x but the full dependence on {x,x') is needed, e.g., to determine the QCD 
evolution of g2 or to compute the high-^r tails of TMDs. Note also that since the 
integrand in Eq. ([TT]) is antisymmetric in x, x', we obtain the non trivial property 







dxg'^{x)^0. (12) 



THE WW RELATION 

Eliminating g"j ^ from ©-dV]) one can derive the Wandzura-Wilczek relation ^ for 
the structure function g2 = —gi + j'LaST' explicitly write down its breaking term 

82{x)-gr{x) = ll^el(tAx)- f^-^gUy)+ [' ^Hy)] , (13) 
^ a V Jx y Jx y J 

Note that g2 explicitly satisfies the Burkhardt-Cottingham sum rule Jq g2{x) = 0, which 
is not in general guaranteed in the OPE [ 12, 13.1. 

A natural question is: how much is the WW relation broken? Model calculations 
have been used to repeatedly argue that the pure twist-3 terms are not necessarily 
small [13|, llTi]. However, in the recent past, since the LIR-breaking gr term was not 
considered in Eq. (fT3l) and the quark-mass term with hi was neglected, the breaking of 
the WW relation was considered to be a direct measurement of the pure twist-3 term 
gT- Therefore, the presumed experimental validity of the WW relation, which we are 
presently going to challenge, was taken as evidence that gr is small. This observation 
was also typically generalized to assume that all pure twist-3 terms are small. 

Our present analysis shows instead that, precisely due to the presence of gr, the 
measurement of the breaking of the WW relation does not offer anymore the possibility 
of measuring a single pure twist-3 matrix element, nor to generically infer its size. On 
the theory side, the quark-target model of Refs. [|l7l. ll8|l can be used to determine both 
gT and gT, which are both comparable in size to the the other twist-2 functions. On the 
experimental side, we used data on polarized DIS on proton and neutron targets to fit the 
WW breaking term A{x) defined as the difference of the experimental data and g^"^' 

Aix)=gf{x,Q^)-gf'^{x,Q^). (14) 

gY"^ was determined using the LSS06 leading twist parametrization [19], and A 
fitted to a functional form allowing for a change in sign and satisfying the Burkhardt- 
Cottingham sum rule. The result is presented in Fig.[Tl and Table[Il where the deviation 
from the WW relation is quantified for a given [x™", jc™^''] interval by 

2 _ /vmin dyxAf^x) 
J/mi„ dyxg^ix) 

with y = log(x). The value of r is a good approximation to the relative magnitude of A 
and g2, which are sign-changing functions. For the proton, we considered three intervals: 
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FIGURE 1. Top panels: the experimental proton and neutron g2 structure function compared to gY^ ■ 
The crosses are computed at the experimental kinematics. The lines are g^^ computed at the average 
of the E155x experiment: the solid (dashed) line is computed with the LSS2006 fits of ^i, with (solid) 



and without (dashed) the HT contribution obtained in the fit. Data points for the proton target 120 
have been slightly shifted in x for clarity. For the neutron only the high precision data from I2lll22 
have been included. Bottom panels: The WW-breaking term Ath for model (1) and (11) compared to the 
higher-twist contribution to . See text for further details. 



TABLE 1. Results of the 1 -parameter fits of the WW breaking term Ath for different 
choices of its functional form, value r of the relative size of the breaking term is computed 
for the whole measured x range, [0.02,1]; the low-x region, [0.02,0.15]; the large-jc 
region, [0.15,1]. See text for further details. 





proton 




X^/d.o.f. not How m 


(I) 
(11) 


Ath =0 

Ath =a{l-x)P{{l5- 
a =0.13 ±0.05 
J3 =4.4±1.0 


h2)x-l) 


1.22 

1.05 15-32% 18-36% 14-31% 


neutron 


(I) 
(11) 


Ath =0 

Ath =a{l-x)P{{l5- 
a =0.64 ±0.92 
J3 =24 ±10 


h2)x-l) 


1.66 

1.11 18-40% 



the whole measured x range, [0.02,1]; the low-x region, [0.02,0.15]; the large-x region, 
[0.15,1]. For the neutron, due to the limited statistical significance of the low-x data, we 
limit ourselves to quoting the value of r for the large-.x; region, [0.15,1]. 

In summary, we have found that the experimental data are compatible with a substan- 
tial breaking of the WW relation in the 15-40% range. 



A PROPOSAL FOR AN EXPERIMENTAL CAMPAIGN 



Figure [T] clearly shows the need for better precision in g2 measurements with both proton 
and neutron targets. In particular, for the neutron high precision is needed away from 
X ~ 0.15 — 0.20 where JLab EOl-012 data almost completely determine the presented 
fits. But even if in the future the WW approximation is found to be more precise than in 
our analysis, we would only be able to conclude that 



This can clearly happen because either gr and gr are both small, or because they 
accidentally cancel each other. Therefore no information can be obtained on the size 
of the twist-3 quark-gluon-quark term gr from the experimental data on g2 aloneH 

However, individually determining the size of gx and gj is very important to gather 
information on the x, x' dependence of the quark-gluon-quark correlator. This can be 
experimentally accomplished by using the EOM ^ and LIR ([7]) and measuring the 

gy^ function, accessible in semi-inclusive deep inelastic scattering with transversely 



In TMD factorization, g\j is the first transverse moment of a twist-2 TMD. Its experi- 
mental determination is challenging because it requires measuring a double spin asym- 
metry in semi-inclusive DIS up to rather large hadron transverse momentum. Further- 
more, in the LIR it appears differentiated in x, which requires a rather fine x binning. 
Preliminary data from the E06-014 and SANE (E-07-003) experiments at Jefferson Lab 
will soon be available, and will demonstrate the feasibility of the proposed measurement 
of and gj. 

This measurement is also very important because the EOM LIR (|7]) and WW 
relation breaking (fT3]) provide 3 independent measurement for 2 independent quanti- 
ties. Verifying them will constitute a pretty stringent test of TMD factorization and its 
connection to coUinear factorization. 
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^ Note that these results were essentially already obtained in Ref. ll24ll . In that work, however, the authors 
assumed gT small and the WW relation small, obtaining a small gj, which is unjustified as we have just 
discussed. 
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